Production of multipartite entanglement for electron spins in quantum dots 



O 
O 
(N 

> 
O 



CD 



o 

> 
O 

o 



X 



F. Bodoky and M. Blaauboer 
Kavli Institute of Nanoscience, Delft University of Technology, Lorentzweg 1, 2628 CJ Delft, The Netherlands 

(Dated: February 1, 2008) 

We propose how to generate genuine multipartite entanglement of electron spin qubits in a chain of 
quantum dots using the naturally available single-qubit rotations and two-qubit Eleisenberg exchange 
interaction in the system. We show that the minimum number of required operations to generate 
entangled states of the GHZ-, cluster and W-type scales linearly with the number of qubits and 
estimate the fidelities of the generated entangled cluster states. As the required single and two- 
qubit operations have recently been realized, our proposed scheme opens the way for experimental 
investigation of multipartite entanglement with electron spin qubits. 
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I. INTRODUCTION 

A. Multipartite entanglement 

Bipartite entanglement refers to non-classical correla- 
tions between two quantum particles, and multipar- 
tite entanglement to non-classical correlations between 
three or more quantum particles. The characterization 
and quantification of the latter is far less understood than 
for bipartite entanglement 0]. In particular, in case of 
multipartite entanglement it is no longer sufficient to ask 
if the qubits are entangled, but one needs to know how 
they are entangled as there are different ways - known 
as entanglement classes - in which three or more qubits 
can be entangled. For three qubits, there are two dif- 
ferent equivalence classes of genuine tripartite entangle- 
ment Q, for four qubits already nine or eight '6], 
and the number of classes is growing with the number of 
qubits. Two entangled states belong to the same equiv- 
alence class and are called locally equivalent if it is pos- 
sible to transform between them using local operations 
and classical communication (LOCC) only, i.e. without 
interactions between two or more qubits. The two classes 
of entanglement for three qubits are the GHZ- and the 
W-class [1, 01 , with representative members [8j] \GHZ^) 
= ^(lOOO)-f |111)) and |W^3> = ^(|001) + 1010) + 1100)) 
(the subscript indicating the number of involved qubits), 
which are both going to be addressed in this paper. Both 
of these classes can be generalized to arbitrary num- 
bers of qubits. Another interesting class of multipar- 
tite entanglement for four or more qubits is the cluster- 
class [9| , which forms the basis of proposals to implement 
a measurement-only quantum computer, the one-way 
quantum computing scheme ^IQ] . These states maximize 
mutual bipartite entanglement and its four-partite rep- 
resentative is 104) = i(|0000)-H 10011) -h|1100)-|llll)). 

Multi-qubit entanglement is thus not a straightforward 
extension of bipartite entanglement and gives rise to new 
phenomena which can be exploited in quantum infor- 
mation and quantum computing processes. For exam- 
ple, there are quantum communication protocols that re- 
quire multi-party entanglement such as universal error 



correction quantum secret sharing ^T^. and tele- 

cloning Also, highly entangled multipartite states 

are needed for efficient quantum computing - all known 
quantum algorithms (such as Shor's factorization [l4i | and 
Grover's search flSsI algorithm) work with multipartite 
entanglement - and GHZ states can be used to construct 
a universal quantum computer [l6j . In addition, multi- 
qubit entangled states provide a stronger test of local 
realism [l^l which is based on individual (rather than 
statistical, as in the bipartite case) measurement results. 
As a general rule, one can say that the more particles 
are entangled, the more clearly non-classical effects are 
exhibited and the more useful the states are for quantum 
applications. 

So far, multipartite entanglement has been realized in 
a number of experiments, using liquid-state NMR , 
photons Jp, H HI, 111, cold atoms Q 111] and 
ions [25l . l26l |. The latter two experiments for trapped ions 
have demonstrated the deterministic creation of a GHZ- 
and a W-state. Tripartite, and more generally multipar- 
tite, entanglement has not yet been realized for qubits 
in a solid-state environment. The latter type of qubit 
systeins, consisting of e.g. electrons confined in quantum 
dots [231 or superconducting Josephson junctions [28j . 
are attractive since they are in principle scalable to an 
arbitrary number of qubits. A number of ideas have 
been suggested for the creation of tripartite entangled 
states, using exciton states in coupled quantum dots 0, 
electron- hole entanglement in the Fermi s ea 13011 and su- 
perconducting charge and flux qubits [sil, l32l. l33l| . 

In this paper, we present schemes for determinis- 
tic creation of GHZ, W- and cluster states for elec- 
tron spin qubits in quantum dots using the naturally 
available two-qubit (Heisenberg exchange) interaction 
and single-spin rotations. This choice of system is 
motivated by the fact that both single-qubit rotations 
and tunable two-qubit Heisenberg exchange interactions 
have already been demonstrated experimentally for these 
qubits [33. 135II. However, our scheme can easily be used 
for other types of qubits as well, e.g. superconducting 
qubits for which tunable coupling has also very recently 
been realized [11]. We show that the required num- 
ber of two-qubit interactions for the generation of N- 
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partite entangled states and for the transformation from 
an disentangled to a maximally entangled basis scales lin- 
early with N for all types of entangled states considered 
here. We also present arguments that the total number of 
single- and two-qubit operations that our schemes predict 
is in fact the minimum number required to create these 
multipartite entangled states using single-qubit rotations 
and Heisenberg exchange interactions. 



B. Electron spin qubits 

An electron spin qubit [3X] consists of a single electron 
confined in a quantum dot (QD), an island in a semicon- 
ducting nanostructure [l^]. The electron occupies dis- 
crete energy levels in the quantum dot which split into 
separate levels for spin-up and spin-down due to Zeeman 
splitting when the quantum dot is placed in an external 
magnetic field. The qubit is encoded in the spin degree 
of freedom, with the ground state spin-up (denoted as 
I t) and defined along the direction of the magnetic field, 
which we assume to be the z-axis) corresponding to the 
logical bit |0) and spin-down (| j)) corresponding to the 
logical bit Electron spin qubits are attractive candi- 
dates for quantum computing since they are in principle 
scalable, relatively robust against decoherence (as com- 
pared to e.g. charge qubits), and allow for a high level of 
control over individual qubits [38^ . 

Two basic operations are available to manipulate the 
state of the qubit: First, coherent rotation of a spin 
around an axis in the {x, ?;)-plane using electron spin res- 
onance (ESR), which consists of applying an oscillating 
time-dependent magnetic field B{t) in this plane whose 
frequency is on resonance with the transition frequency 
between | "f) and | |). A rotation around a certain angle 
is controlled by the time of application and the strength 
of the magnetic field, and is described by the evolution 
operator 



and we will denote it from now on as 



UR(t) = exp {ij/2) / B{t) k-adr 



(1) 



which corresponds to the Hamiltonian 7i_R,(t) — 
— {l/2)h'yB{t)k- a, where 7 denotes the gyromagnetic ra- 
tio, a = {ax, ayyCTz), and k = (sin 9 cos (f>, sin 9 sin0, cos 9) 
represents a unit vector on the Bloch sphere [0e[O, tt), 
(f>£[0, 2tt)] in the direction of the magnetic field. The evo- 
lution operator ^ corresponds to a rotation i?^"-* (/3) = 
exp(— (l/2)j/3fc • f?^"^) of the n-th qubit with angle /3 = 
—7 /p B{t) dr around axis fc, where tr^"^ = 1 (g) . . . (g) ct (g) 
. . . 1. These ESR-induced rotations have recently been 
experimentally observed in quantum dots (33 | . 

The second available operation is Heisenberg interac- 
tion between two spins described by the evolution oper- 
ator Usxit) = exp{-ih/4) /o J(t) ct^") • dr. Here, 
J(r) is the time-dependent exchange energy. By tuning 
the interaction time t with a gate voltage, the (SWAP)"- 
gate can be directly generated as (SWAP)" = UEx{t), 



(^^SWAP)'^ 





cos(^7r) 
i sin(^7r) 






i sin(^7r) 
cos(^7r) 










(2) 

where a{t) = —j J{T)dT. For a — ^ the y/UswAP gate 
maximally entangles two spins of opposite directions. A 
vC/swAP-operation has also recently been demonstrated 
for spin qubits 35]. 

Together, single-qubit rotations and the V^^SWAP-gate 
form a universal set of quantum gates, into which any 
quantum operation can be decomposed [37| . 



C. Outline 

This paper is organized as follows: in Sec.|IT]we show how 
to generate A^-partite entangled cluster states fSec. lH A]) . 
GHZ-states (Sec HTB]) and W-states (Sec. lHCl) using the 
smallest number of two-qubit (C/swAp)"-operations and 
single-qubit rotations. In Sec. IHII we analyze the effects 
of errors in the timing of the (C/swap)"- and single-qubit 
operations on the generation of A^-partite cluster states, 
quantified by the fidelity. Finally, in Sec. IIVI we discuss 
the feasibility of the multipartite entanglement genera- 
tion scheme that we propose in the context of present-day 
experimental techniques, followed by conclusions. 



II. GENERATION OF MULTIPARTITE 
ENTANGLED STATES 

In this section we describe the generation of A^- 
qubit entangled states in a chain of quantum dots 
where each dot is occupied by one electron, using 
single-qubit rotations and pairwise exchange interac- 
tions between nearest-neighbor spins. Starting from the 
ground state which consists of N disentangled up-spins 
|00...0) and using a recursive approach in A^ we de- 
rive sequences of single-qubit rotations and two-qubit 
(C^swAp)"-operations which, when appfied to |00...0), 
yield a A^-partite cluster-, GHZ- or W-state. We be- 
gin by briefly recounting the generation of entanglement 
and the implementation of a basis transformation for two 
qubits, and then present our main results in Sees. Ill Al 
inU and HIC] below. 

For two qubits, the shortest sequence required to trans- 
form the ground state |00) into a maximally entangled 
state is 



(3) 



with k an arbitrary axis in the (a;-y)-plane and z = 1, 2. 

It has also been shown that the shortest sequence re- 
quired to implement the transformation from the stan- 
dard (or computational) basis to a maximally entangled 
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basis consisting of Bell states is given by [36 



_g(l,2,i) _ 
k 



(1,2) 
SWAP 



(1,2) 
SWAP 



(4) 



with i — 1,2. Two VC^swAP-operations are needed in 
since one C/s wap -interaction only entangles two of the 
four standard basis states. Given a linear array of quan- 
tum dots in which each dot is occupied by a single spin 
qubit in the |0) or |1) state, entangled states of three or 
more qubits can be generated by pairwise application of 
the sequence ([4]), as we show in the next subsection. By 
applying local operations in between, one can control to 
which class the generated entangled states belong. This 
forms the basis of our calculations in the next two sub- 
sections. Without loss of generality, we choose the axis 
k as X and i = 1 in Eq. and omit the indices k and i 

in E^j^'^'^"^ in the following. 



A. Cluster states 

It is straightforward to see that pairwise application of 
Eq. (j4|) to a chain of N disentangled qubits, each of which 
are either in the |0) or |1) state, results in an A^-partitc 
entangled cluster state, as we prove in Proposition III. II 
below. Specifically, we will prove that application of the 
sequence 



transforms an arbitrary disentangled state of the iV- 
partite standard basis into a cluster state. The definition 
of linear cluster states for iV qubits is as follows: the clus- 
ter state is the state resulting when applying the Ising in- 
teraction Z("'"+i)(6') = exp(-i6'/4) (l-ai"^)(l-cri"+^^) 
with 6 — TT (the so-called z-phase gate) to each 
neighbor in a A^-qubit chain prepared in the 
state (g)^il/\/2(|0) + The z-phase gate 

can be generated in quantum dots as = 



{n,n+l) 
SWAP 



(see [131), where we omit the overall phase factor (as in 
the rest of this paper). In order to generate a cluster 
state, one thus has to apply this z-phase gate to each 
pair of qubits in the state 1/V2(|0) + |1)) = i?;;(7r/2)|0). 
With these observations, we are now ready to prove 

Proposition II. 1. The two sequences 

^ z^N^l,N) ^(1,2) (g^) 

and 



(6b) 



Proof. We start by rewriting the sequence as: 

ZN ^ z(^-l,W)i?W(^)...z(2.3)i?(3)(^)Z(l-2) 

"2 "2 
xi?f(f)4'n|) (7a) 

= (7b) 

^ 2 

which consists of iV — 1 applications of the operator 
^(«,n+i) ^ z("^"+i)i?.^"+^^(7r/2), plus an additional ro- 
tation Riy\'K/2). We now write in terms of 
^(n,«+i) jgq^ ^-^^Yi the goal to express as = 
LE^i, where L is a product of local operations. To this 
end, we use the identity 



Z 



(1:2) 

(8) 

and substitute Eq. ([5]) into Eq. (j7bp . The two rotations 
on qubit 1 commute with the sequence to the left of them, 
and can thus be absorbed into the local operation L, leav- 
ing Z(2.3) acting on R^y\TT/2)RP {~tt/2). By rewriting 

z("-"+i)i?(")(|)4")(-|) 

= (9) 



repetitive substitution of Eq. (|9|) into Eq. (j7b|l for in- 
creasing n and using commutation relations to reorder 
the resulting sequence such that rotations are shifted to 



(5) the left of all \/?7swAP-operations, wc find 



with 



Z^ = LE^^-'^^K..E^''^^ =LE^i, 



(10) 



(11) 

Since L does not change the entanglement class of the 
state that has been generated by E^^, we have thus 
proven that application of the transformation E'^ to any 
state of the standard basis leads to a cluster state. □ 

Inspecting the sequence (j6bp . we see that 2{N — 1) 
VC/swAP-operations and {N — 1) rotations are needed to 
generate a A^-partite cluster state, whereas when using 
previously proposed implementations of the z-phase gate 
^(«,n+i) ^ ^ ^Q^g^j q£ _ 3^ rotations are required. 
The sequences given by Eqs. (jS]) transform the standard 
basis into a basis of cluster states. Note that in order to 
transform the ground state |0 . . . 0) into a cluster state, 

the sequence . . . E^^''^') ^fu^J^Ri^\Tr), which 



are locally equivalent. 



contains one v^^swap operation less than (j6bp . is suffi- 
cient, since \/u^^lp Rx'\t^) [see Eq. already maxi- 
mally entangles the first two qubits. 
A special property of the z-phase gate, being a diagonal 
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matrix, is that the matrices commute for dif- 

ferent n. As a result, the transformation ([5a|) from the 
standard basis to the cluster basis can be done in two 
steps: first, all the even-numbered qubits are simulta- 
neously entangled to their (odd) neighbor to the right, 
and then the same is done for the odd-numbered qubits. 
In the next proposition, we show that the same commu- 
tation relation applies for the (non-diagonal) 
matrices in Eq. ([5]). 

Proposition II. 2. The two sequences 

^^(N-2,N-1) ^(2,3)^W(^) ^(1)(^) (^2a) 

and 



^(N-l,N)^{N~?.M-2) ^ ^ ^ ^(1,2) 
^ ^(iV-2.iV-l) _g(2,3) 



(12b) 



are locally equivalent. 

Proof. Since the 2''^"'"+^) matrices commute for differ- 
ent n, the sequence (|12ap is equivalent to the right-hand 
side of Eq. (|6ap and hence by Proposition lII.il to the se- 
quence (j6bp . It now remains to be shown that (j6bp is 
equivalent to (|12bp . This directly follows from the fact 
that the matrices and i;("+i.«+2) [Eq. g])] com- 

mute: 



= 0. 



(13) 



The two sequences (|12al) and (I12bp differ by the same 
overall local operation L [Eq. lfTT]) ] as in Prop. UlT] □ 



Propositions III. II and III. 21 imply that the sequences 
(j6bp and (jl2bp . which consist only of the entangling 
operations create a state of the cluster class. 

Can this be done with less operations? To answer this 
question, consider first three qubits: we have already 
seen that the first -\/C/swAP-operation in the sequence 
([6b)) can be omitted when entangling the ground state, 
and if we start in an appropriate excited state also the 
first rotation is not needed. That leaves the sequence 

f^swAP-R^^^ \J t^swAP \l t^swAP ■ It is Straightforward to 
check, e.g. by calculating the tangle r [12] of the re- 
sulting entangled state, that this is the shortest sequence 
of VC^SWAP-operations and single qubit rotations that 
creates a tripartite cluster state: if any of the four op- 
erations is omitted, r = and the resulting state is no 
longer a cluster state. Generalizing to an arbitrary num- 
ber of qubits, we note that omitting any operation in the 
sequence (j6bp or (jl2bp leads to a state which is not max- 
imally connected in the sense of Ref. Q, and therefore 
cannot be a cluster-state. 



B. GHZ states 

In this section we show how the disentangled N- 
qubit state |00...) can be transformed into a A^-qubit 



GHZ-state using single-qubit rotations and yC/swAP- 
operations. We start with the observation that GHZ- 
states are generated by successive application of the 
CNOT-gate C/cnot: 

Observation 1. Starting from the disentangled N- qubit 
state I00...0), the N-partite GHZ-state \GHZn) = 
|00 • • • 0) -|- |11 • • • 1) (disregarding normalization) is gen- 
erated by N — 2 applications of the C/cnot • 



\GHZn) = 



n 



u, 



{n,n+l) 
CNOT 



xi?i''(|) A^wap4'H^)|00---0)'(14) 

where C^cnot^'' denotes a Uqnot- operation with the n-th 
qubit as control bit, and the {n-\- l)-th qubit as target bit. 

Note that the order in the product in Eq. ITU) , starting 
with the highest n = iV — 1, is essential. In Eq. the 



operation4^)(-f )i?ii)(f )i?i'^(f ) V^C/^^wap W on 
the first two qubits yields the Bell state (l/^/2)(|00) -\- 
|11)) and each successive t/cNOT-gate entangles one more 
qubit to this superposition, resulting in the A'^-partite 
GHZ-state \GHZn). Using single-spin rotations and 
VC^SWAP-operations only, the shortest sequence of op- 
erations required to implement the t/cNOT-gate is given 
by El 



r(l,2) 



CNOT 



(15) 



We now substitute (fTS]) into Eq. (HH). By moving all 
single-qubit rotations that commute with the sequence 
of operations to the left of them in front of all vSWAP- 
operations and defining 

(16) 



Eq. (|T4l) becomes 



GHZn) = L 



n 



U, 



(n,n+l) 
CNOT 



'^/u^Ri'\n)\00■■■0) N>3, 



(17) 



where L consists of single-qubit rotations. We see from 
Eq. pT]) that a total of (2A^ — 3) \/?7swAP-operations 
and a minimum of (3A'^ — 5) single-qubit rotations are 
needed to transform the separable state |00 . . .) into an 
A^-partite entangled state in the GHZ-class. Compared 
to cluster states (see the previous section) , the generation 
of a GHZ-state thus requires {2N — 4) more single-qubit 
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rotations. In practice, implementation of the sequence 
P?)) can be done in the most efficient way by applying two 
-y/CTswAP-operations simultaneously in each step. This 
can be achieved by starting with qubit number m = N /2 
(for N even, or m = {N + l)/2 for N odd) in the middle 
of the chain and reordering the sequence ([TT]) as (for N 
even) : 

m+l , — — ^ 

nl TT(N-j+l,N-j) jA],j+l) \ 
I'^CNOT '^CNOT ) 

]=N-l ^ ^ 

x^t^i^^^'WlOO-'-O)- (18) 

The two C/cNOT-operations in between the brackets in 
(jlSp can be performed simultaneously. An analogous ex- 
pression as (dHJ apphes if N is odd. 
Using the same line of reasoning as in the previous sub- 
section and the fact that GHZ-states are also maximally 
connected [;9j, one can directly show that the number of 
VC^swAP-operations and the number of rotations over tt 
in (|17p is minimal. Although we have no formal proof 
for this, we suspect that the total number of rotations in 
P7p is minimal. 

C. W-states 

A W-state for N qubits is a multipartite entangled state 
that can be written in the form 

\Wn) = ^(|10 . . . 0) + |01 . . . 0) + . . . + |00 . . . 1)). 
yN ^ ' 

TV terms 

(19) 

In this section we show how \ Wm) can be generated using 
the least possible number of single-qubit rotations and 
VC^swAP-operations when starting from N disentangled 
qubits in the state |00 . . . 0). Since the W-state (|19p con- 
sists of a superposition of 7V-qubit states in which one bit 
value differs from all the others (representing for exam- 
ple a collective spin state of N spins with one distributed 
excitation), a straightforward "recipe" to generate a A^- 
partite W-state starting from the ground state |00 . . . 0) 
is to flip one qubit to the value " 1" and then distribute 
this bit value equally over all qubits such that each qubit 
is excited with a fraction 1/N . Rotating qubit 1 and then 
applying a (?7swap)" interaction to qubits 1 and 2 results 
in the fraction of the first excitation being cos (a/2); so 
in order to have the first qubit excited with a fraction 
1/A'', a has to be chosen as a = 2 arccos(-y/l/A^). By 
applying the same reasoning to all the following qubits, 
we find 

\Wn) = 



with /i„ ~ 



where (^) denotes a (?7swAp)''-gate on the qubits 

n and n -|- 1 . 

For implementation of the sequence (|20p we can apply 
the same trick as for the GHZ-states by starting from 
the middle qubit in the two opposite directions along 
the chain and performing two VC^SWAP-operations si- 
multaneously. To find the correct interaction angles we 
then need to distinguish between even and odd number 
of qubits N. For the even case, we rotate the qubit 
m = N/2, apply an interaction J^'^-^+i) (1/2) and then 
proceed as with two independent strings of length m. For 
the odd case, we rotate the qubit m = (A^+1) /2, apply an 
interaction j(™''"+^)(/z) with /i = 2 a,rccos{ y/m/N) and 
then proceed as with two independent strings of length m 
(to the left) and m — 1 (to the right). Thus the required 
number of operations to transform the state |00...0) 
into a A^-partite W-state consist of (A^ — 1) (?7swap)^- 
operations and one rotation. 



III. FIDELITY 

In the previous sections we have assumed perfect control 
of the single- and two-qubit operations, i.e. we assumed 
that all the pulses were perfectly timed. However, for 
a physical implementation it is important to estimate 
the effect of imperfections in the (C/swap)" gate oper- 
ations and in single-qubit rotations on the intended fi- 
nal entangled states. In this section, we provide such 
an estimate for the generation of cluster states [Eq. ([5])], 
assuming that the control of each VC^SWAP-operation is 
off by a small parameter e, and similarly for each rota- 
tion by a small parameter S, i.e. we replace VC^SWAP 
(t^swAp)^5+0 and R{tt) R{tt + 5), where |e| < 1/2 
and <C tt. As measure for the effect of the inaccura- 
cies e and 6 we use the fidelity F [i^, which describes 
the overlap between the intended ("perfect") state \(j)N) 
and the generated (in the presence of the inaccuracies) 
state \4>n)- For pure states, F is defined as: 



F=^{j0Ar|0A,)(0]v|0Jv), (21) 

where |0Ar) and \4>n) are both normalized. 

In the following, we calculate F as a function of A^ for 

the cluster state 

\cj)r,)^E^''-^'''K..E(-^'^^\00...), (22) 



j(Af-i,Af)^^^_^-j • • • J^^' i?^^^(7r)|00- • •0)iip to second order in e and S. We start by expanding 

(t/swAp)^+' [Eq. (El)] and R^in+S) [Eq. (P)] up to second 
(20^rder in 6 and e: 
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(C^SWAP) 





i+e_^-^(i+e) COs(f(i- 

jsin(f(i 




)) 
^)) 



t sm 
cos I 





.2^2 ' 





^)) 







-(3+«^) 



(23) 



Rx[7T + S] = 









cos(^) 


-zsin(^) 



-I sm [-^^Y^ 


cos(^) 




Since we assume |e| and \S\ to be small compared to 1/2 and tt, respectively, we can expand Eqns. 
second order in e and S: 



and 



(24) 



to 



ITT ZTT 



cos 



(4 + 2- 



ITT \ / ZTT 

exp [ — 1 exp I ye 



l + iA-^ + 0(A3) 



. /TT TT \ . /7r\ /TT \ / TT \ . / TT 

sm 1 e] = sm — cos — e + cos — sm — 1 

14 2/ 14/ V2/ V4/ V2 



1 

71 



1 

V2 



TT 

2 



cos I — ) COS I — e 



sm — sm — e 



l-A- 



= sm 



TT 

2 

52 



cos 



0(^3) 



COS 



sm 



l-^ + 0(i?3), 



= -B + 0{B^ 



(25a) 



(25b) 



(25c) 



(25d) 



Here c = exp(i7r/4) is a constant, A = 7re/2, and i? = 6/2. Using Eqs. 
operation E = (f/swAp)^^*^ Rx'\'^ + 5) (C/swap)^^"^ up to second order in the parameters e and 5: 



(25e) 

5]), we construct the entangling 



E = 



-c(l + 2iA)B ^ 
^[l + (l + z)^-(l-z)^2_^] 

^[1-(1-,:)A-(1+,:)A2_^] 



;)5[l + (l + z)A-(l-z)A2 
-2c^AB 



-cB 







V2 



|[l-(l-z)A- (1 + ^)^2 



2 



Si 

2 



^[l-(l-z)A-(l + z)A2 
-cB 
~2c^AB 
^[l + (l + z)A-(l-z)A2 



The cluster state \4>n) is defined as 



(27) 



2 



Si 
2 





^[l~(l-.)A-(l+.)^2_^] 
^[l + (l+.)yl- (1-^)^12-4!] 

-c(l + 2iA)B 



(26) 



r 



The order of the entangling operations in 
Eq. (|27p has to be the same as for the intended state 
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FIG. 1: (a) 3D-plot of the fidelity F, Eq. (f28)l . for iV = 3 as a 
function of the two parameters e and 5. The red line and the 
transparent red plane show the cut along which the 2D-plots 
in part (b) are taken, (b) The decrease of the fidelity F when 
the number of qubits A'^ is increased. We chose here 5 — e/2. 



|0Ar) in Eq. (f^^ . since different ordering generates dif- 
ferent states (although in the same entanglement class). 
We now calculate \(t)N) and \(j)N) and from these the fi- 
delity F [Eq. (|21l)] for an increasing number of qubits. 
We then find Hjl 



F = \ 1- {N - 1)^2 



5iV-9 



> 3, (28) 



for (A^- 1)^2 - (5iV- 9)^2/2 < 1. We have verified that 
(PS)) is valid up to = 10, and suspect that it is true for 
all values of A^. Remark that the fidelity F increases as 
the square root of N, the number of qubits. Figure ((T^) 
shows the fidelity as a function of the inaccuracies e and S 
for the case of three qubits. Numerical evaluation shows 
remarkably high fidelities even for systems with many 
qubits, e.g. F{N = 10,e = 0.05, S = 0.1) = 0.95. 
This suggests that the proposed sequences for the cluster 
states enable generation of many-qubit entangled states 
with high fidelity. Although the fidelity of generating 
e.g. A^-partite GHZ-states (not shown) are lower because 
there are more single-qubit rotations required to trans- 
form the A^-qubit ground state into a GHZ-state, they 
are also within reach of experimental implementation, as 
we discuss in the next section. 



IV. DISCUSSION 



spins IS 

quires 




tron spins in quantum dots. As demonstrated experimen- 
tally, the duration of a VC^swAP-operation of two electron 
35| . and a spin rotation over tt/2 re- 
As a rough estimate, we then find 
that the sequence Eq. (jl2bp to implement a Af-partite 
cluster state using simultaneous application of y/UswAP- 
operations to all pairs of qubits requires ~ 2-50 = 100 ns 
for any A^ (since the linear array of qubits can be entan- 
gled in just two steps: first each of the even-numbered 
qubits to their right neighbor, and then the same for the 
odd-numbered qubits). The time required to implement 
the A^-partite GHZ state (fTT]) depends on the number 
of qubits and amounts to ^ • 100 ns (using the 

ordering given in Eq. (jlSp ). The limiting time for the 
implementation of these sequences of operations is the 
decoherence time T2, which has not yet been measured 
for a single spin. Rabi oscillations of a single electron 
spin [s^l have been seen for more than 1 /is, indicating a 
decoherence time T2 > 1 ^s. Based on this estimate for 
T2, the generation of a A^-partite cluster state and a A^- 
partite W-state thus seems feasible for any A^ in a time 
shorter than T2, whereas the generation of GHZ states 
should be possible for up to ~ 10 qubits p^ . 
To conclude, we have calculated general sequences to 
generate genuine A'^-partite entangled states in various 
entanglement classes starting from a separable A'^-qubit 
state in the computational basis and using the least pos- 
sible number of single-qubit rotations and two-qubit ex- 
change interactions. For all entangled states that we con- 
sidered (cluster states, GHZ-states and W-states) we find 
that the total number of operations required to gener- 
ate these sates scales linearly with the number of qubits 
N. The generation of A^-partite W-states requires the 
least amount of operations, namely (A'^ — 1) vC/swAP- 
operations and 1 rotation. They are followed by the A^- 
partite cluster states that require a minimum of {2N ~ 3) 
exchange interactions and (A^ — 1) rotations and the GHZ 
states that also require a minimum of (2A'^ — 3) exchange 
interactions and {3N — 5) single-qubit rotations. We also 
calculated the fidelity F for the generation of A^-partite 
cluster states in the presence of imperfect single-qubit 
rotations and Vt^swAP-operations, and find that F de- 
creases as F ~ ~ M-^ + as the number of qubits 
grows, with fi,i> > and fiN — v <1. 
Our results can be implemented for electron spins in 
quantum dots [4^, for which Heisenberg exchange is the 
naturally available two-qubit interaction in the system. 
We estimate that our proposed scheme for the generation 
of multipartite entangled states is feasible for at least 10 
qubits within current experimental accuracy. Finally, we 
emphasize that the approach used in this paper can be 
used for any kind of two-qubit entangling interaction and 
provides an analytical scheme to calculate the implemen- 
tation of multipartite entangled states for any type of 
qubit. 



In this section we briefly discuss the experimental feasi- 
bility of generating multipartite entangled states of elec- 
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